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Abstract 

The phase shift of the 0(4) symmetric </> 4 theory in the symmetric 
phase is calculated numerically using the relation between phase shift 
and energy levels of two-particle states recently derived by Liischer. 
The results agree with the prediction of perturbation theory. A prac- 
tical difficulty of the method for a reliable extraction of the phase shift 
for large momenta due to the necessity of a precise determination of 
excited two-particle energy levels is pointed out. 
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Monte Carlo simulation techniques enable us to study various non-perturb- 
ative phenomena of field theories defined on a space-time lattice. Because of 
the Euclidean formalism of lattice field theories, however, the method does 
not allow a direct evaluation of physical quantities characterizing scattering 
processes such as phase shift . 

It has been known for some timed, |, that finite -size effects in energy 
levels are closely related to scattering amplitudes. Indeed, asymptotic be- 
havior of energy levels of two-particle states for large volume can be written 
in terms of scattering length || (see Refs. ||, ||, |6], [7] for its applications). 
Recently this asymptotic formula has been generalized to an exact relation 
between energy levels of two-particle states in a finite box and phase shift 
[§]. Since two-particle energy levels are calculable through standard Monte 
Carlo techniques, this relation opened a possibility of extracting phase shift 
through numerical simulations. Of particular interest is that the relation 
may allow a determination of resonance parameters in QCD[||. 

The corresponding relation in two dimension has been used to numeri- 
cally extract the phase shift of the 0(3) non-linear a model[10| and a coupled 
Ising system [|TT| and a good agreement has been found between the numer- 
ical results and analytic predictions. These results encourage us to examine 
the practical applicability of the method for realistic field theories in four di- 
mensions. In this article we report on our attempt to extract the phase shift 
of the four-dimensional 0(4) symmetric 4 theory in the symmetric phase. 
This model provides a good testing ground of the method since the triviality 
of the theory, which is confirmed in many waysfll^] though no exact proof 
exists, enables us to check the results against perturbative calculations. 

Consider a system of two identical particles in a cubic box of a size L 3 , 
whose states are classified by the irreducible representations of the cubic 
group S0(3,Z). Let Wj (j = 0, 1, 2, • • • ; W < W x < ■ - •) be the energy 
levels of the two-particle states in an irreducible representation of SO(3,Z). 
Let m be the mass of the particle and define the momentum kj corresponding 
to Wj by 



Wj = 2^/m 2 + kj. (1) 

Liischer's formula|| relates kj to the set of phase shifts 5i(kj) having the 
angular momentum I not excluded by the symmetry of the two-particle state. 
For the states in the Af representation the relevant phase shifts are Sq, 64, 
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8q, ■ ■ '. At low energies the s-wave phase shift do dominates. Neglecting the 
phase shifts Si with / > 4 the formula takes the form|J 

=^(^)-jV, (2) 

where <p(q) is given by 

Zoo{l;q 2 ) -tTT^ 2 q 
with Z 00 (l;q 2 ) defined by an analytic continuation of 

2 00 ( S;9 2 ) = -L^— i-. (4) 
v% n t^3 (n 2 - q 2 ) s 

In a numerical simulation of a given size L, the formula (|2|) determines the 
phase shift only at a discrete set of momenta k , k±, ■ ■ ■ corresponding to the 
energy levels Wo, W±, • • •. The momenta, however, can be shifted by using a 
different lattice size. Combining the results for various lattice sizes one can 
obtain the full momentum dependence of the phase shift. 

We apply the method above to the 4 theory in four dimensions with 
0(4) symmetry, employing the standard lattice action given by 

s = -2«EE On +A , E = i. (5) 

nfj, a =l a=l 

The system undergoes a second-order phase transition at k = k c = 0.30411(6), 
above which the 0(4) symmetry is spontaneously broken. The simulations 
are made at k — 0.297 in the symmetric phase on an L 3 x 16 lattice with 
L = 8 and 12. We carried out 4 x 10 6 sweeps of the heat bath algorithm for 
each lattice size, measuring observables at every 10 sweeps. 

We extract single-particle energies from the exponential decay of the prop- 
agator 

E<WW)> — const.e-^)* (6) 

a 

with 0p(£) the projection of 0° to the spatial momentum p at a time slice 
t. The exponential fits are made over the range t = 3 — 6 for the spatial 
size L = 8 and t = 4 — 6 for L — 12. The results for the momenta p = 
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(0, 0, 0), (1, 0, 0), (1, 1, 0) and (1, 1, 1) (in units of 2vr/L) are given in Table 1. 
The errors are estimated by the jackknife method with the bin size of 10 6 
sweeps. Our values for the mass m — E(0) are consistent with those of Ref. 
H for the same choice of k and L. 

The 0(4)-scalar two-particle operator for a center-of-mass momentum p 
is given by 

OpW = EEW%W ( 7 ) 

a R 

where the summation over cubic rotation R ensures the projection onto the 
yl^-sector. The corresponding two-point function is defined by 

GVW = (0p(*)0p'(O)> - (O p (t))(O p ,(0)) (8) 

Inserting the complete set of states, one can rewrite the two-point function 

as 

G PP '(t) = E«'e-^, (9) 

j 

with v p the coupling of the state j to the two-particle operator O p . The 
two-particle energy levels W/s can be extracted by diagonalizing G pp >(t) as 
a matrix in p, p' at each time slice t, and fitting the eigenvalues Xj(t) to a 
single exponential |10| 

XM^C-e-W (10) 

In our calculation we truncate the momenta to the subset p = (0,0,0), (1,0,0), 
(1,1,0), (1,1,1). For the lowest state a y 2 fit is made to the exponential form 
( |10|) for t = 2 — 7, while for higher states we used the range t = 2 — 4 as 
errors in the eigenvalues become quite significant for larger values of t. The 
results for the two-particle energies are given in Table 2 where we also list 
the values of the momentum k = \J (W/2) 2 — m 2 with m = E(0) taken from 
Table 1. The errors are estimated by the jackknife method with the bin size 
of 10 6 sweeps. 

We have made a one-loop perturbative calculation of the phase shift to 
compare with numerical results. For the 0(iV)-symmetric 4 theory in the 
continuum defined by the Lagrangian 

1 N 1 N 1 / N \ 2 

C =nJ2 9^ a dy a - -ml £ - 7[«R E ^ + (counterterms), 

* a=l 1 a=l 4 - \a=l I 

(11) 
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the two-particle scattering amplitude up to one-loop level, renormalized by 
the momentum subtraction at p = 0, is given by 



4tt 2 

T(k, k') = —A 



or + l<&{A<p(s) + B{ V {t) + cp(u)) + 0(a R )} 



(12) 



where s, t, u are the Mandelstam variables constructed from the initial and 
final momentum k and k' in the center-of-mass frame, k = Ikl = Ik'l, ckr, = 



^r/167t 2 , andcufc = \Jk 2 + m R . The coefficients take the values A = (iV+2)/3 
and B = 1 for the scalar channel (^4 = 2/3 and B = (A^ + 6)/6 for the tensor 
channel), and the function (p is defined by 
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(p(z) = / dxln 1 T x(l - x) . (13) 

Jo V m R J 

The phase shifts Si (k) are defined through the partial- wave expansion of 
T(k,k') 

Stt 00 p 2i5 ; (fe) _ i 

T(k,k') = £(2l + l)fl(cosg) , (14) 

with # the angle between k and k'. The s-wave phase shift 5 (k) is therefore 
obtained as 

= ~— A[a K + a 2 R {Af(k) + Bg(k)} + 0(4)], (15) 
z to k 

where the functions f(k) and g(k) are given by 

„/, x k , LUh + k , 

m = — in - — i, i6 

Uk m R 

,(*) = ^sm^^ln 1 

+ ir^sm^^4ln^^ + fcCOSi -l). (17) 
2 Jo y k cos § m R ' 

Liischer's formula @) is derived in the continuum space-time and the 
dispersion relation of the one-particle energy E(p) = Vp 2 + rn 2 enters into 
the proof in an essential way. The simulation results for the one-particle 
energy are compared with the continuum dispersion relation in fig.l. We 
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find a good agreement (E(p) / a/p 2 + m' 2 =0.994(6)) for the lowest momentum 
p = 2tt/L ~ 0.5 for the lattice size L = 12. An increasing deviation for larger 
momenta indicates that finite lattice spacing effects become non-negligible 
for higher momenta. 

In fig. 2 we show by filled circles our results for the s-wave phase shift 5o(k) 
in the scalar channel extracted from the lowest two-particle state (see Table 
2 for numerical values). The dotted lines show the function <ft(kL/27r) and 
the solid line represents the one-loop perturbative result fllS]) calculated with 
the infinite volume estimates for the renormalized parameters = 0.3044 
and «r = 0.142 for k = 0.2970. We have also used the data of ref.|| 
for the lowest two-particle energy to calculate the phase shift. The results 
are plotted by open circles in fig. 2 and are consistent with our values. We 
observe a reasonable agreement between the results of simulations and that 
of one-loop perturbation theory. A trend may be present, however, that 
the numerical results become smaller than the perturbative prediction as the 
lattice size decreases. This may represent a systematic bias due to an increase 
of vacuum polarization effects for small lattice sizes, which are not taken into 
account in (0). 

As is seen from the figure, the phase shift is almost linear in the region 
of momenta which can be explored by the lowest two-particle states. These 
states therefore do not give information on the phase shift beyond the scat- 
tering length, which has already been calculated d, |5j using the leading term 
of Liischer's formula in 1/L expansion derived earlier [Q. In order to extract 
the phase shift for large momenta we must therefore examine excited two- 
particle states. The rightmost point in fig. 3 shows the phase shift extracted 
from the first excited state for k = 0.297 and L = 12. The agreement with 
the perturbative prediction, albeit with a sizable error, is quite encouraging. 

We note that fig. 3 reveals an important point in practical applications 
of Liischer's formula. The function (f)(kL/2ir) increases very rapidly for mo- 
menta corresponding to excited states (this is specific to four dimensions; 
in two dimensions 0(A;L/27r) = kL/2 and the slope is constant ||10||). As a 
consequence, a reliable extraction of the phase shift for large momenta re- 
quires quite a precise determination of the momentum k = <J (W/2) 2 — m 2 , 
and hence that of the two-particle energy W. For example, reducing the 
error of the phase shift for the first excited state in fig. 3 to a 10% level ne- 
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cessitates a calculation of the momentum within 0.5% accuracy compared 
to the 2% error of the present data . Achieving such an accuracy requires 
an order of magnitude more computer time (our run for the size L — 12 
took 20 hours for 4 x 10 6 sweeps on HITAC S820/80 with the peak speed 
of 3 Gflops). The requirement becomes rapidly more demanding as we go 
to higher excited states. For instance, with the present statistics we only 
obtained — 5o(k) = 0.4(9) for the second excited state. 

To conclude, the method works well for momenta which can be explored 
by the lowest two-particle state. We find, however, that extraction of phase 
shifts for large momenta requires quite a high statistics determination of 
excited two-particle state energies due to the steep increase of the function 
(j)(kL / 2-n) . We feel that this presents an obstacle in practical applications of 
the method, with which one can in principle extract scattering data beyond 
scattering lengths in four dimensions. 

I would like to thank Professor A.Ukawa for leading me to this interest- 
ing subject and for continuous support. I am also thankful to Professors 
M.Fukugita, M.Okawa and Dr.H.Mino for kind advice. The computations 
were carried out on HITAC S820/80 at KEK. 
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P (0,0,0) 


(1,0,0) 


(1,1,0) 


(1,1,1) 


L = 8 0.321(2) 
L = 12 0.305(2) 


0.819(7) 
0.602(3) 


1.11(2) 
0.783(6) 


1.24(4) 
0.91(1) 



Table 1: Single-particle energy E(p) in lattice unit for k = 0.297 for L = 8 
and 12. 



L J 





1 


2 


3 


8 W j 


0.73(2) 


1.6(2) 






kj 


0.17(2) 


0.74(11) 






-Sofa) 


0.14(4) 








12 ^ 


0.646(9) 


1.26(3) 


1.66(9) 


1.72(15) 


kj 


0.106(6) 


0.552(11) 


0.77(7) 


0.80(10) 


-Sofa) 


0.12(2) 


0.40(15) 


0.4(9) 





Table 2: Two-particle energy Wj, the corresponding momentum kj in lattice 
unit and the s-wave phase shift — So(kj) calculated through (2) for k = 0.297 
for the size L = 8 and 12. Values for the cases when errors are too large are 
omitted. 
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Figure captions 



Fig. 1 -E(p)/Vp 2 + rn 2 as a function of p — |p| for L = 8(squares) and 
L— 12 (circles) at n = 0.297 with m = 0.321(2) for L = 8 and m = 0.305(2) 
for L = 12 as input. 

Fig. 2 S-wave phase shift calculated from the lowest two-particle state en- 
ergy for k = 0.297 with various lattice sizes L. Momentum k is in lattice 
unit. Dashed lines represent the right hand side of (2) for each L. Solid line 
represents the one-loop result (15). 

Fig. 3 S-wave phase shift calculated from the first excited two-particle state 
energy for k = 0.297 and L = 12. Those from the lowest ones are also plotted. 
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